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Many  studies  have  investigated  the  role  of  near-surface  bubbles  in  ambient  noise 
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generation,  acoustic  backscatter,  and  transmission  loss.  Medwin,  Novarini  and  Bruno, 

3 

and  Hall  have  addressed  the  impact  of  bubbles  for  the  special  case  of  surface  duct 
propagation.  Hall  has  developed  a  comprehensive  semi-empirical  acoustic  model  that 

4  5 

incorporates  the  bubble  measurements  of  Johnson  and  Cooke  and  Thorpe.  Hail's 
model  lends  itself  toward  studying  the  impact  of  a  bubbly  layer  on  the  general  surface 
loss  problem.  Transmission  loss  for  acoustic  environments  that  involve  many  surface 
interactions  (e.g.,  shallow  water  and  surface  duct)  is  particularly  affected  by  propaga¬ 
tion  loss  at  the  near-surface  boundary.  Thus,  the  inclusion  of  near-surface  bubbles  in 
the  acoustic  propagation  problem  is  seen  as  a  necessary  requirement  in  ensuring  that  all 
mechanisms  at  the  surface  boundary  are  properly  addressed.  This  study  will  examine 
the  impact  of  the  bubbly  layer  on  acoustic  propagation  in  the  0.5-  to  5-kHz  frequency 
regime. 
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BACKGROUND 


•  A  HOST  OF  STUDIES  HAVE  SHOWN  THE  INFLUENCE  OF 
NEAR-SURFACE  BUBBLES  ON  AMBIENT  NOISE,  ACOUSTIC 
BACKSCATTER  AND  TRANSMISSION  LOSS.  HALL  (J. 
ACOUST.  SOC.  AM.  VOL  86,  NO.  3,  SEPTEMBER  1989)  HAS 
DEVELOPED  A  COMPREHENSIVE  ACOUSTIC  MODEL  BASED 
ON  EMPIRICAL  MEASUREMENTS  OF  BUBBLE  POPULATIONS 
INVESTIGATING  ONE-WAY  SURFACE  DUCT  PROPAGATION 
LOSS  AS  A  FUNCTION  OF  SOURCE  DEPTH. 


•  HALL'S  MODEL  IS  EXTENDED  TO  EXAMINE  THE  INFLUENCE 
OF  THE  BUBBLY  LAYER  ON  NEAR-SURFACE  PROPAGATION 
AND  RESULTING  IMPACT  ON  SURFACE  LOSS  MEASUREMENTS 
AND  MODELING 


2 


Theoretical  studies  of  the  acoustic  reflection  loss  at  the  air-sea  boundary  generally 
neglect  the  influence  of  the  near-surface  bubbly  layer.  The  present  analysis  addresses 
the  impact  of  refraction  and  absorption  on  the  acoustic  wave  just  before  and  after  the 
interaction  with  the  sea  surface.  The  objective  is  to  examine  the  relative  contribution 
of  the  bubbly  layer  to  the  total  near-surface  loss  as  a  function  of  grazing  angle,  wind 
speed,  and  frequency. 

This  study  will  apply  ray  theory  to  the  sound  propagation  near  the  surface.  In  the 
viewgraph,  the  bubbly  layer  is  modeled  to  result  in  (1)  a  depth-dependent  sound  speed 
(C(z))  yielding  a  change  in  grazing  angle  (from  6  to  e  ')  and  (2)  absorption  along  the 
acoustic  path  length. 
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NEAR  -  SURFACE  PROPAGATION  IN  A  BUBBLY  LAYER 
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To  simplify  the  problem,  this  study  is  based  on  the  following  assumptions: 

1.  The  existence  of  a  horizontally  stratified  layer  of  bubbles  without  inhomo¬ 
geneities  in  the  bubble  population  density.  Future  work  in  this  area  will  address 
Langmuir  "banding"  and  bubble  plume  distributions. 

2.  The  absence  of  scattering  effects  due  to  the  frequencies  examined  being  low 

enough  to  result  in  a  very  small  product  for  the  acoustic  wavenumber  k  and 

a 

the  bubble  radius  a. 

3.  No  contribution  to  the  surface-scattered  field  from  off-specular  reflections. 

4.  The  exponential  decay  of  the  sound  speed  anomaly  at  the  surface.  It  will  be 
shown  that  this  assumption  is  supported  by  measurements.  It  has  been  found 
that  invoking  the  first  assumption  results  in  negligible  bubble  effects  below 
10  m.  When  bubble  plume  dynamics  are  considered,  it  is  expected  that  depths 
greater  than  10  m  will  need  to  be  addressed. 

5.  No  orbital  motion  and  other  turbulence  due  to  surface  waves  in  the  bubbly  layer. 
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ASSUMPTIONS 

1.  BUBBLE  DENSITY  IS  AZIMUTH  ALLY  HOMOGENEOUS  FOR  A  GIVEN 
DEPTH.  LANGMUIR  MECHANISM  (BANDING)  IS  NEGLECTED. 


2.  SCATTERING  DUE  TO  BUBBLES  IS  NEGLECTED  (k*  a  «1). 


3.  INCOHERENT,  OFF-SPECULAR  SEA-SURFACE  SCATTERING  IS 
NEGLECTED.  SURFACE  LOSS  MODELS  ARE  ROUGH-SURFACE 
ACOUSTIC  SCATTERING  MODELS. 


4.  DEPTH  DEPENDENCE  OF  THE  SOUND  SPEED  ANOMALY  AND 
ATTENUATION  DUE  TO  BUBBLES  CAN  BE  APPROXIMATED  BY  AN 
EXPONENTIAL  DECAY.  BUBBLE  EFFECTS  ARE  NEGLECTED  AT 
DEPTHS  > 10m. 


5.  THE  EFFECTS  OF  SEA-SURFACE  MOTION  ON  THE  BUBBLY  LAYER  ARE 
NEGLECTED. 
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The  Hall  model  is  extended  to  obtain  closed-form  expressions  for  the  sound  speed 
anomaly  and  attenuation  as  functions  of  frequency,  depth,  and  wind  speed.  The  sound 
speed  anomaly  is  used  with  Snell's  law  to  determine  the  total  refractive  effect,  and  the 
new  acoustic  wave  direction  (0*  )  is  used  in  the  rough  surface  scatter  model  to  be 
implemented.  Rough  surface  scattering  is  treated  as  independent  of  the  bubbly  layer. 
The  total  attenuation  is  found  through  numerical  integration  over  the  acoustic  path 
length  before  and  after  the  surface  reflection.  Small  uncertainties  in  the  near-surface 
boundary  loss  can  result  in  large  transmission  loss  uncertainties  for  multiple  surface- 
interaction  propagation  environments.  The  effects  of  the  bubbly  layer  may  be  small  at 
lower  frequencies  (~0. 5-1.0  kHz),  but  should  still  be  examined  for  accurate  modeling  of 
such  environments. 
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APPROACH 


•  EXTEND  HALL  MODEL  TO  OBTAIN  SOUND  SPEED  ANOMALY  AND 
ATTENUATION  AS  FUNCTIONS  OF  FREQUENCY,  DEPTH  AND  WIND 
SPEED.  DERIVE  CLOSED-FORM  EXPRESSIONS. 

•  NUMERICALLY  INTEGRATE  OVER  ACOUSTIC  PATH  LENGTH  TO 
DETERMINE  TOTAL  REFRACTION  AND  ATTENUATION  AS  A  FUNCTION 
OF  GRAZING  ANGLE  AND  WIND  SPEED. 

•  DETERMINE  "EFFECTIVE”  SURFACE  LOSS  DUE  TO  BUBBLE 
ATTENUATION  AND  ROUGH  SURFACE  SCATTERING. 

•  PERFORM  MODEL  VALIDATION  OF  SURFACE  LOSS  AND 
PROPAGATION  LOSS  WITH  EMPIRICAL  DATA  COLLECTED  IN  SURFACE 
DUCT  AND  SHALLOW  WATER  ENVIRONMENTS. 
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The  complex  sound  speed  in  the  bubbly  layer  (C)  is  found  by  evaluation  of  the 

integral  for  all  bubble  radii  (a).  Hall  models  the  bubble  population  density  spectrum 

4  5 

level  by  an  empirical  fit  to  the  Johnson  and  Cooke  measurements  and  Thorpe.  The 
real  part  of  C,  Re(C),  yields  the  sound  speed  anomaly,  and  the  imaginary  part  of  C, 
Im(C),  is  used  to  compute  the  attenuation.  Damping  occurs  primarily  from  thermal 
diffusion  for  the  low  frequencies  of  this  study.  ^ 
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COMPLEX  SOUND  SPEED 


(D 


where 


N(a)  =  BUBBLE  POPULATION  DENSITY  SPECTRUM  LEVEL 
Co  =  SPEED  OF  SOUND  WITHOUT  BUBBLES 
6  =  DAMPING  COEFFICIENT 
a  r  =  RESONANT  RADIUS  (AT  FREQUENCY  f) 
a  =  BUBBLE  RADIUS 
f  =  FREQUENCY  (kHz) 
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The  sound  speed  anomaly  (Re[C]  -  CQ)  is  expressed  in  equation  (2).  The  n/T  and  w 

dependence  is  obtained  directly  from  reference  3.  Hall  noted  that  the  depth 

dependence  of  the  anomaly  is  approximately  an  exponential  decay.  The  depth 

dependence  is  found  to  have  a  decay  constant,  m(w),  that  is  dependent  on  wind  speed. 

The  wind  speed  dependence  shown  in  equation  (3)  was  determined  using  a  least-squares 

fit  to  the  sound  speed  anomaly  generated  from  equation  (1)  for  the  range  of  wind  speeds 

(0-20  m/s)  at  1.25  kHz.  A  check  of  the  validity  of  the  below  expression  is  done  by 

7 

comparison  with  the  measurements  of  Farmer  and  Vagle. 
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SOUND  SPEED  ANOMALY 

HALL  INTEGRATION  RESULTS  SHOW  THAT  DEPTH  DEPENDENCE  OF 
THE  SOUND  SPEED  ANOMALY  CAN  BE  APPROXIMATED  BY  AN 
EXPONENTIAL  DECAY: 


Re[C]  -C0  = 


where 


m(w)  =  .0078w2  -  0.33  w  +  4.36 

w  =  WIND  SPEED  (m/s) 
z  =  DEPTH  (m) 


(2) 


(3) 


COMPARISON  WITH  FARMER  AND  VAGLE  (JASA,  86(5),  NOV  1989) 
AT  5  kHz  SHOWS  REASONABLE  AGREEMENT. 
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Measurements  by  Farmer  and  Vagle  are  done  at  5  kHz  for  two  different  wind  speed 
conditions  (10  m/s  and  12-15  m/s).  Reasonable  agreement  between  the  da*a  and  the 
approximation  to  Hall  is  seen  for  both  wind  speeds.  The  sound  speed  anomaly  data  at 
the  lower  wind  speed  are  shown  to  split  to  — 5  m/s  and  -3  m/s  at  the  surface  (z  =  0). 

The  data  collected  by  Farmer  and  Vagle  are  represented  by  the  lower  values  ( — 5  m/s). 


and  the  exponential  fit  to  the  data  derived  by  reference  7  results  in  the  higher  values 
(-3  m/s). 
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The  rate  of  attenuation,  a,  is  found  from  equation  (4)  after  evaluation  of  the 
complex  sound  speed.  The  approximation  to  a  yields  a  quadratic  dependence  of 
frequency  and  an  exponential  dependence  of  depth  with  a  decay  constant  of 
(l/2)('s/w/15)(l/L(w))l  where  L(w)  is  the  e-folding  depth  of  the  bubble  population  density 
spectrum  level  N(a).  L(w)  is  empirically  fit  by  Hall  and  given  by  the  bilinear  curve 

L(w)  =  0.4  (m)  w  <  7.5  m/s 

0.4  +  0.115  (w  -  7.5)  w>  7.5  m/s. 

The  frequency  dependence  differs  slightly  from  Hall  due  to  the  inclusion  of  lower 
frequencies  (0-5  kHz)  in  the  quadratic  approximation  (least- squares  fit). 
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RATE  OF  ATTENUATION 

a  =  -(ln?o)2nf  He) 

(4) 

OBTAIN  a  FROM  HALL  INTEGRA!  iON  RESULTS.  FIT  FREQUENCY 
DEPENDENCE  WITH  QUADRATIC  POLYNOMIAL  AND  APPROXIMATE 
DEPTH  DEPENDENCE  WITH  EXPONENTIAL  DECAY: 

a  =  (.0061  f2  +  .2066  f  -  .0521  )  (-^-)3  e(2  V 

,W5  L(w)j  •  ‘5> 

WHERE  a  IS  IN  dB/m. 
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The  rate  of  attenuation  showing  the  exponential  depth  dependence  is  plotted  at  a 
frequency  of  1.25  kHz.  It  is  noted  that  the  rate  of  attenuation  is  relatively  small  for 
even  large  wind  speeds,  but  when  the  total  attenuation  is  found  by  integration  over  the 
acoustic  path  length,  the  impact  can  be  significant  at  low  grazing  angles. 
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Exponential  Depth  Dependence  of  Attenuation 
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The  total  or  "effective"  near-surface  loss  is  modeled  to  consist  of  the  following 
two  independent  components: 

1.  The  rough  surface  scattering  term  (SL  ),  where  the  grazing  angle  utilized  (e'  ) 
is  a  result  of  refractive  effects  of  the  bubbly  layer  and 

2.  The  integrated  attenuation  loss  term  (SL  ),  where  the  rate  of  attenuation,  a,  is 

EL 

integrated  over  the  acoustic  path  length  before  and  after  the  interaction  with 
the  surface. 
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"EFFECTIVE”  NEAR-SURFACE  LOSS 


SL  =SLr  +SLa 


m 


where  SLr  =  ROUGH  SURFACE  REFLECTION  LOSS  WITH 
INCIDENT  GRAZING  ANGLE  e'g 

INCORPORATING  REFRACTION  EFFECTS 
DUE  TO  BUBBLES. 


6'g  =COS-1 


COS 


(eg) 


fRe(C)  -Cp 


+  1 


SLa  =  INTEGRATED  ATTENUATION  LOSS  DUE  TO  A 
BUBBLY  LAYER. 

r- ° 

=  2-  a(z)csc[e(z)]dz 

Jz/C0 
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The  wind  speed  dependence  of  refraction  is  shown  by  comparison  of  the  grazing 

angle  of  the  acoustic  wave  at  the  bottom  of  the  bubbly  layer  (0  )  with  the  grazing  angle 

at  the  surface  (©'  ).  At  very  low  grazing  angles  associated  with  surface  duct 
B 

propagation  (e  <  2°),  it  is  seen  that  even  a  moderate  wind  speed  of  10  m/s  can  triple 
8 

the  grazing  angle  to  ~6°.  Shallow  water  propagation  can  have  grazing  angles  that  are 
typically  less  than  10°,  and  it  is  seen  that  at  moderate  wind  speeds  (15-20  m/s),  there  is 
a  significant  increase  in  the  grazing  angle  relative  to  the  bottom  of  the  bubbly  layer. 
This  plot  shows  that  even  in  the  absence  of  attenuation,  refractive  effects  alone  will 
result  in  a  flattening  of  the  grazing  angle  dependence  of  the  near-surface  reflection 
loss.  Similar  results  are  seen  in  reference  8,  where  the  grazing  angle  distribution  of  an 
acoustic  pressure  field  shifts  to  higher  angles  when  a  bubbly  layer  is  introduced. 


VIEWGRAPH  12 
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A  power  law  fit  is  applied  to  the  integrated  attenuation.  The  result  provides  a 
closed-form  expression  that  brings  out  the  dependence  on  wind  speed,  grazing  angle, 
and  frequency.  The  limit  of  application  of  the  equation  is  for  ©  >  2°  and  w  <  20  m/s. 
The  quadratic  fit  for  f  does  not  exhibit  appreciable  error  above  -0.5  kHz. 
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INTEGRATED  ATTENUATION  (SL  ) 
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INTEGRATE  ATTENUATION  OVER  ACOUSTIC  PATH  LENGTH  BEFORE 
AND  AFTER  SURFACE  REFLECTION  AT  AIR-SEA  BOUNDARY.  RESULTS 

CAN  BE  APPROXIMATED  BY  A  POWER  LAW  FIT  IN  W  AND  0. 


SLa  =(.0375  f2  +  1 .2968  f  -.320)  •  0a  ■  10b  +  c 

where  a  =  .0019  w2  +  .1348  w  - 1.8612 
b  =  -.0047  w2  +  .1869  w- 1.24 
c  =(.0009  w2  -.0483  w  +  .3582)  (log(e))2 

6  =  grazing  angle  (deg) 
f  =  frequency  (kHz) 
w  =  wind  speed  (m/s) 
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The  integrated  attenuation  due  to  bubbles  (SL  )  is  plotted  for  both  the  numerical 

a 

solution  and  the  power  law  fit.  The  interesting  result  is  the  increase  in  attenuation 
with  decreasing  grazing  angle  and  increasing  wind  speed.  Although  the  integrated 
attenuation  is  small  for  the  1.25-kHz  example  shown,  the  loss  can  be  seen  to  be 
significant  for  many  surface  interactions  associated  with  a  surface  duct  or  shallow 
water  environment.  Also,  relative  to  rough  surface  scattering  theories  (e.g.,  Kirchhoff. 
perturbation  .  .  . ),  where  the  surface  reflection  loss  approaches  zero  at  zero  degrees 
grazing  angle,  the  SL  is  an  important  term  to  reconcile  in  modeling  and  measuring 
near-surface  propagation. 


NUWC 


Integrated  Attenuation  Due  to  Bubbles 


A  comparison  of  the  integrated  bubble  attenuation  (SL  )  is  made  with  the  loss 

<1 

associated  with  the  limiting  ray  cycle  of  the  surface  duct  measurements  of  Marsh  and 
9  10 

Schulkin  as  implemented  by  Weinberg.  The  data  and  the  SL&  modeling  results  show 
that  at  moderate  to  low  wind  speeds  (w  <  ~8  m/s),  the  Marsh-Schulkin  loss  can  be 
explained  exclusively  by  the  mechanism  of  attenuation  in  a  near-surface  bubbly  layer. 
There  is  an  ~3-dB  difference  between  SL&  and  the  empirical  data  for  w  >  10  m/s.  This 
difference  suggests  that  the  modeling  needs  to  address  either  (1)  a  loss  at  the  air-sea 
boundary  (SLp  at  higher  wind  speeds  or  (2)  a  greater  attenuation  loss  (perhaps  due  to 
bubble  clouds  or  plumes). The  attenuation  loss  implies  that  the  stratified  layer 
assumption  of  the  model  may  need  to  be  relaxed  at  the  higher  wind  speeds. 


VIEWGRAPH  15 


16 


12 

The  modified  Eckart  surface  loss  model  was  selected  as  an  example  of  including 
SLf  in  the  predictions  of  the  total  "effective"  near-surface  loss.  The  Marsh-Schulkin  data 
are  plotted  from  0°  <  6  <2°  for  wind  speeds  of  5,  7.5,  and  10  m/s.  The  "effective"  loss, 
SL,  is  plotted  for  e  >  3°  both  with  and  without  bubbly  layer  effects.  A  flattening  of  the 
grazing  angle  dependence  at  lower  grazing  angles  is  seen,  and  the  modeled  SL  approaches 
the  Marsh-Schulkin  empirical  results  for  wind  speeds  below  -10  m/s.  Thus,  the 
discrepancies  between  the  low  grazing  angle  empirical  data  and  the  rough  surface 
scattering  theory  (below  10  m/s)  can  be  explained  by  inclusion  of  the  bubbly  layer 
refraction  and  attenuation  mechanisms  in  the  surface  interaction  problem. 
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The  acoustic  model  of  Hall  that  addresses  the  effect  of  bubbles  on  near-surface 
propagation  has  been  extended  to  a  study  of  the  generalized  surface  loss  problem.  The 
modeled  sound  speed  anomaly  in  the  bubbly  layer  is  shown  to  agree  with  measurements 
at  5  kHz  by  Farmer  and  Vagle.  The  integrated  attenuation  due  to  the  bubbly  layer  is 
shown  to  increase  with  decreasing  grazing  angle  and  increasing  wind  speed.  The 
refraction  effects  result  in  an  increase  in  the  grazing  angle  at  the  surface  that  becomes 
more  pronounced  at  higher  wind  speeds.  The  "effective"  near-surface  loss  at  low 
grazing  angles  will  be  greater  than  sea-surface  interface  scattering  without  bubbles. 
The  "effective"  near-surface  loss  is  shown  to  be  similar  to  the  Marsh-Schulkin  surface 
duct  data  at  1.25  kHz  for  moderate  to  low  wind  speeds  (<10  m/s). 
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SUMMARY 

1.  MODELING  RESULTS  OF  SOUND  SPEED  ANOMALY  IN  A  BUBBLY  LAYER 
SHOW  AGREEMENT  WITH  MEASUREMENTS  AT  5  k HZ  BY  FARMER  AND 
VAGLE.  THE  IMPACT  OF  THE  SOUND  SPEED  ANOMALY  IS  TO  INCREASE 
THE  INCIDENT  GRAZING  ANGLE  AT  THE  SURFACE  RELATIVE  TO 
BUBBLE-FREE  WATER. 


2.  THE  INTEGRATED  BUBBLE  ATTENUATION  INCREASES  WITH 
DECREASING  GRAZING  ANGLE  AND  INCREASING  WIND  SPEED. 
ATTENUATION  HAS  A  QUADRATIC  DEPENDENCE  ON  FREQUENCY  WITH 
NOTABLE  EFFECTS  ABOVE  0.5  kHZ. 


3.  ATTENUATION  AND  REFRACTION  EFFECTS  DUE  TO  A  BUBBLY  LAYER 
RESULT  IN  AN  INCREASE  IN  THE  EFFECTIVE  NEAR-SURFACE  LOSS  AT 
LOW  GRAZING  ANGLES  RELATIVE  TO  SEA-SURFACE  INTERFACE 
SCATTERING  WITHOUT  BUBBLES.  LOW  WIND  SPEED  RESULTS  ARE 
SIMILAR  TO  MARSH-SCHULKIN  MEASUREMENTS  AT  1.25  kHZ. 
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